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ABSTRACT 

A differential method of orbit improvement utilizing observational 
data is presented for Vinti^s spheroidal solution of the dynamical 
problem of unretarded artificial satellite motion about an oblate planet, 
recently modified so as to permit the exact inclusion of the effects of 
the third zonal harmonic term of the planet’s gravitational potential 
field. The first-order Taylor’s series expansion used for the equations 
of condition is fitted to observational values by an' iterated least- squares 
process, producing successively smaller corrections to the orbital 
elements. A mean set of elements, conditioned by the observations, re- 
sults for use in orbital predictions for intermediate time points or for 
later epochs. The differential coefficients in the conditional equations, 
applicable to any type of observational data, are derived analytically 
from the equations of the accurate reference orbit. The method of dif- 
ferential correction may be used for orbits of all inclinations, including 
equatorial and polar; both a first-order treatment and a lengthier, but 
more exact, second-order treatment for the periodic variables are 
given. 
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DIFFERENTUL CORRECTIONS APPLIED TO VINTI’S 
ACCURATE REFERENCE SATELLITE ORBIT WITH INCLUSION 
OF THE THIRD ZONAL HARMONIC 

by 

Harvey Walden and Stan Watson 
Goddard Space Flight Center 


INTRODUCTION 

In order to predict the precise position, at a given time, of a satellite revolving about a planet, 
an extensive mathematical theory of the satellite's motion and the exact values of certain physical 
parameters (e.g., gravitational constants) are required. If the position of an artificial satellite of 
the earth is to be determined relative to an observer on the earth’s surface, then, in addition, the 
accurate geodetic position of the observer is necessary. Any mathematical theory of motion is 
based upon certain constants of the motion which initially must be determined empirically. Once 
approximate values for these constants of the motion are available, then they may be utilized in the 
theory to predict future orbital positions of the satellite. A comparison between the positions pre- 
dicted by the theory and those actually observed will indicate discrepancies, expressed numerically 
as the differences between the observed positions and the respective computed positions. These 
differences are known as observational residuals, and the fact that they are nonzero is due to 
several causes. One fundamental cause which unfortunately is ever present is the inadequacy of 
the theory, in that it is unable to account mathematically for all of the physical forces acting upon 
the satellite. Moreover, there are always errors associated with the observations themselves, be- 
cause of fluctuations in the atmospheric density, optical imperfections in the telescope, inaccura- 
cies in the reduction of the observational recordings, and the like. Further, the physical parameters 
required by the theory often are not known with sufficient precision, and the geodetic position of the 
observer likewise may be measured inaccurately. Finally, there are errors introduced by the 
approximate values for the constants of the motion used in computing the theoretical position. 

Despite all of the inadequacies and shortcomings enumerated, the very existence of the dis- 
crepancies between theory and observation provides a means of improving the approximate values 
of the constants of the motion. One procedure for improvii^ these constants in an analytic theory 
of motion involves expressing the incremental changes in the positional coordinates due to changes 
in the constants of the motion as coefficients having the form of partial derivatives in a Taylor’s 
series expansion. In such a case, the partial derivatives must be determined as e^licit functions 
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of variables which arise in the analytic theory of motion. This method of orbit improvement 
is known as differential correction. 

The purpose of this paper is to provide such an orbit improvement method for the spheroidal 
theory of artificial satellite motion. The spheroidal method for satellite orbits, developed by Vinti 
(References 1 through 5), supplies a procedure for calculating an accurate reference orbit of any 
dra^-free satellite moving in the gravitational field of an axially symmetrical oblate planet. In 
the case of artificial satellites of the earth, the intermediary reference orbit reproduces exactly 
the zeroth and second zonal harmonic coefficients in the series expansion of the geopotential func- 
tion, and also accounts for more than half of the fourth zonal harmonic. Recently, Vinti (Refer- 
ences 6 and 7) has modified the spheroidal potential to allow exact inclusion of the effects of the 
third zonal harmonic as well, heretofore the major neglected effect in the spheroidal orbital theory. 
(The first zonal harmonic is eliminated entirely, of course, by proper choice of coordinate origin.) 
Accounting for the third zonal harmonic directly in the intermediary orbit in this manner affords a 
more accurate treatment (Reference 8) than would be possible through perturbation theory. This 
paper will present a differential method of orbit improvement based upon the modified spheroidal 
theory that includes the effects of the third zonal harmonic term. 


FUNDAMENTAL EQUATIONS FOR DIFFERENTIAL CORRECTION 

Consider a system of constants of the motion q. (i = 1, 2, . . . , n) which are utilized in a math- 
ematical theory of motion to predict orbital positions of a satellite. In this context, such constants 
of the motion are generally referred to as orbital elements, and the number n contained in the 
system is often 6. Denote a positional coordinate of the satellite at time t by R( t ) , using sub- 
scripts and to distii^ish between observed values and values computed by the analytical 
theory using the orbital elements. We assume that the differences of the coordinates, R^Ct) - R^(t) , 
as well as the required corrections to the elements in the improvement of the orbit, are sufficiently 
small so that their squares and higher powers may be neglected. We can then express the obser- 
vational residuals by a truncated Taylor^s series expansion restricted to first powers in the general 
form 


, (?R^ (t) 

R,(t)-R^(t) - ^ -■ Aq. . 

i=i 

Here the computed value R^(t) is viewed as a function of n independent variables, q.(i = 1, 2,>- ,n), 
which are to be improved by the additive increments Aq. (i = 1, 2, • • • , n ). The coefficients in the 
Taylor’s e3q>ansion e>q)ressing the increment of the coordinates caused by a change in the orbital 
elements have the form of partial derivatives, and these must be determined anal 5 rtically from the 
equations in the mathematical theory of motion. Each separately observed coordinate yields an 
equation for corrections of the elements of the form given earlier. Ordinarily, the available number 
of such so-called equations of condition far exceeds the number n of unknowns Aq. . This set of 
linear simultaneous equations forms an inconsistent system (because of inherent random, and 
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possibly systematic, errors in the observations) for which no exact solution exists. It may be 
solved by the method of least squares, yielding ^’preferred^’ values for the unknowns. The solution 
of the equations by the principle of least squares follows well-known schemes (Reference 9), so 
that the problem reduces to evaluatii^ the derivatives of the coordinates with respect to the 
elements (the so-called differential coefficients). 

We now introduce certain conventions in order to allow discussion of the problem in more 
explicit terms. Assume that the satellite's positional data are recorded at the tracking stations in 
the form of direction cosines observed with respect to a topocentric (i.e., situated at the earth*s 
surface) coordinate system. The topocentric or ’local*' coordinates will be distinguished by the 
subscript, ”m”. The system is orthogonal and right-handed, with the plane tangent to the 

earth’s surface. The Xj,-axis extends in an easterly direction along the line of latitude; the Y^,-axis 
extends in a northerly direction along the line of longitude, and the z„-axis is normal to the surface 
and points toward the geodetic zenith. If and denote the observed direction cosines of a 
satellite (in the and Y^-directions, respectively) for a given time of observation, then the 
correspondii^ computed values of the direction cosines are given in terms of the local coordinates 

by 


Xm 

The computed value of the third direction cosine, (in the Zj^,-direction), is not an independent 
parameter, but is predetermined by and through the relationship 

For this reason, each satellite observation provides two, and only two, independent coordinate 
values, which here are chosen to be L and M . 

The mathematical theory of motion will ordinarily predict a satellite’s position with respect 
to a rectangular geocentric system (x, Y , z ) of coordinates, the so-called inertial frame of refer- 
ence. In this system, the X - Y plane is the earth’s equatorial plane, and the origin is situated at 
the earth’s center of mass. The X-axis extends toward the vernal equinox (the first point of Aries); 
the Y-axis is orthogonally to the east to form a right-handed system, and the z-axis coincides with 
the earth’s polar axis. In order to obtain a satellite’s local coordinates from its inertial coordi- 
nates, the inertial coordinates of the observation point at the time of observation, which we shall 
denote (x^, y^, Z^) , must be known and two rotations performed to bring the local and inertial 
systems into parallel alignment. Let us denote by the ai^le between the vernal equinox and the 
X„-axis as measured in the observation latitude plane. Then will depend upon the loi^itude of 
the tracking station and the hour angle of the vernal equinox at the time of observation. If repre- 
sents the latitude of the trackii^ station, then the local coordinates of the satellite at observation 
time are given by the matrix equation 
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1 

0 

0 " 


cos 

sin 

0 



x" 




0 

sin 0JJ 

cos 


- sini/f^ 

cos 

0 



Y 

- 

Yt 


0 

- cos 0 ^ 

sin 0 ^ 


0 

0 

1 



z 


Zt 


— - 


— 




— 



; — 


— 



The difference of column matrices on the extreme right represents a translation from the earth's 
center to the tracking station position; the center matrix on the right represents a rotation in the 
latitude plane about the polar axis through an angle to brii^ the inertial x-axis into coincidence 
with the station's X^^-axis; the remaining matrix on the right represents a rotation in the longitude 
plane about the -axis through an angle equal to the complement of 0 ^^ to brii^ the inertial z-axis 
into coincidence with the station's -axis. This matrix equation, when expanded, reads 

X„ = (X-X^cosv^^ + (y-Y^ sin^, , 

Yj, = - (X ~ Xj. j sin sin + (y - Y^ ) cos sin + (Z - Zj ) cos , 


and 


~ ( X - X^) sin cos - ( Y ~ ) cos cos + (Z - Z^) sin . 

We can now write the first-order Taylor's series expansion for the equations of condition 
corresponding to each time of observation in the following more explicit form: 

Lo - = 

and 


Mo - = 

Since the local coordinates are functions of the inertial coordinates, which are in turn functions 
of the orbital elements, then the chain rule may be used to expand the differential coefficients. Thus 



5Lc 



■?Y„ 

3K 9Z„ 

aq. 

■ 

dq- 


dqi 

^ <?Zm 


and 











aqj 


5q^ 


<9qj 



i=l 
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The partial derivatives of the direction cosines with respect to the local coordinates are found 
directly from the expressions for andM^. Thus 


and 


dX^ 






XmYm 


i9Z„ 


XmZm 




XmYm 


aK 

av„ 




V 


aM^ 


YmZm 

V 


where 


Rm - ( Xm ^+ Y „^. Z „^)>/2 


Since the coordinates x^, Y^, and and^the angles 
ments (and merely geodesic functions), we have the matrix equation 


and 0^ are independent of the orbital ele- 


^X„ 


aqj 


5Y„ 


aqj 







sin 


- cos 


jin (9- 


cos 0^ 


- sin0 


sin 0^ 


cos 0^ 


0 


0 


dX 


0 


dY 


dZ 


Equivalently, 


^Xm _ ax . , aY 

^ V/, + sin./-, -5^ > 
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dX 

- sin./-^sin0Dai: 




and 

, r, ^ 5Y . „ dZ 

aqT - sin.//^ cos " c°s cos ^ • 

The problem remains to evaluate the differential coefficients , aY/dqj , and a^aqj (i = 1, 2, 

• • • , 6), which are the derivatives of the inertial coordinates with respect to the orbital elements. 
At this point, the mathematical theory of motion becomes of primary importance. 

METHOD USED TO EVALUATE DIFFERENTIAL COEFFICIENTS 

The constants of the motion q. (i = 1, 2, • . . , 6), which we choose for the mean orbital 
elements for the modified spheroidal theory of satellite motion (including the exact effects of the 
third zonal harmonic coefficient of the oblate planet's gravitational field) are the following: 

qj = a, the semimajor axis; 

= e, the eccentricity; 

q^ = S, corresponding to sin^ I in Keplerian motion, where I is the inclination of the orbital 
plane to the equator. (However, s may be negative for orbits sufficiently close to equatorial. For 
a polar orbit, S is unity, and in all cases, < 1.); 

q^ = /Sj , corresponding to the negative of the time of passage through pericenter in Keplerian 
motion; 

qj = /Sj , corresponding to the argument of pericenter in Keplerian motion; and 

qg = /Sj , corresponding to the ri^t ascension of the ascending node in Keplerian motion. 

These elements differ slightly from those selected by Vinti in the final algorithm for the reference 
orbit (Reference 7, Section 12), but the foregoing parameters seem more suitable for the differen- 
tial correction. 

The relationship between rectai^lar inertial coordinates and oblate spheroidal coordinates 
P, V) and 4>, which are involved in the solution of the problem of satellite motion, is given by 

X = y(p* + C* ) ( 1 - Tj^) cos <i> , 
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Y 


j/(p^+c^)(l-77^)sin0 , 


and 


Z - p7) - S , (-1 ^T) <l) . 

Here c and S are adjustable parameters that are chosen to agree with the coefficients of the zonal 
harmonics in the series eiq)ansion of the earth’s potential function. In terms of the earth’s equa- 
torial radius and the zonal harmonic coefficient the proper choices are 

For the earth, the values are c = 210 km and 8=1 km. 

In order to evaluate the derivatives of the inertial coordinates with respect to the orbital 
elements, we must know dp/dq^y and d<P/dq^ (i = 1, 2, • - • , 6 ). The oblate spheroidal co- 

ordinates are rather involved functions of the orbital elements. The process for determining the 
partial derivatives of p, 77 , and is a lengthy one which will be presented in a synthetic, rather 
than analytic, manner. That is, necessary partial derivatives of the simpler functions of the orbital 
elements will be given first, followed by partial derivatives of more complicated functions of the 
elements involving the predetermined partial derivatives. 

It is worth noting in passing that the equations to be presented herein apply to orbits of all 
inclinations, including equatorial. There are no special simplifications introduced in the case of 
equatorial or near-equatorial inclinations, as there were for the spheroidal theory that did not 
include the effects of the third zonal harmonic (Reference 5 ). 

However, as in the earlier spheroidal theory, the differential correction may include derivatives 
of the periodic terms taken through the second order, or, alternatively, it may be simplified to 
omit periodic terms higher than first order. In either case, some second-order effects are in- 
cluded with the first-order terms and even with the zeroth-order terms. The terms of second and 
higher orders that are to be omitted in the simplified version will be indicated as such. 


TIME-INDEPENDENT PARTIAL DERIVATIVES IN THE DIFFERENTIAL CDRRECTIDN 

We begin with Equation 23.3,* which reads 


p - a(l -e*) . 


•Until otherwise indicated, all equation numbers used in specifying the defining relation for a given variable refer to Reference 7. 
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Thus, 


da 


^ = 1 - 


i£ 

de 


= - 2 ae 


and 


<5P _ . 

■5B' - 0 . 

Since = -A/2 , by line 8 of Section 12, Reference 7, then by Equation 40.1 we find 

ac^ (ap - S) (1-S) - 4a^ S(1 ~S) [l + c^ (ap)"^ (3S - 2)] 

^ (ap - c^) (ap-c^s) + 4a^c^ S + 4c^ S^p“^ S(l-S) (3ap-4a^ - c^) 

If we denote the numerator and denominator of bj by N and D respectively, then we have 


(9bj 

da 


= - r2 


(l-S)D-i |c2s-a(2p+a|f) + 4S^s[| (2 - f|f ) + (f )' (3S - 2) (l - 2 f )_ 


- ND"^ <( 2 ap-c^-c^S)^p + a^^ + 8 ac^ S 


+ 4 ( 1 ) S2S(1-S) 3p-8a+ ^8 +2 ^ -3aj ll 


5bj 

i^e 


(1-S) fill + 4(1)' s[] 


1 2 ^ (3S - 2 )J I' - ND'2 |£ a(2 ap - S ) 


/C\2 / o2 ^2 

+ 4(p) S2S(1-S) (8 ^ + 2 - 3 


and 


<9bj 

~dS 


2 n-i 


ac^ D 


ap + (1 - 2S) (4 ^ 8^ +c/ -4(^) S2 (gS^ - lOS + 2 ) 


- ND~^ (4 a^ - ap + c^) + 4^p^ 8^(1- 2S) (3 ap - 4 a^ - 


Since by line 8 of Section 12, Reference 7, then by Equation 40.2 we find 
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- bja“^ (ap-c^) , 


Thus 



and 


dh^ 


- \ (®t>2)' 


> db, 
"5S 


This last sequence of three equations can be written as a single generalized equation if we intro- 
duce the following notation. Let = a, q 2 = e, and q3 = S, and let be the Kronecker delta, 
defined as follows: = 1 when i = j, and S.. = 0 whenever i / j . Then, for i = 1, 2, 3, 


(9b2 


1 

2 



- bi (pSii + a 



This notation will be used extensively in order to write later sequences of equations in an efficient 
manner. 

Combining Equations 35.1 and 35.2 with equations already mentioned gives 


Thus 


and 


3Po r ^^2 ( dp\ ( ^^t \] 

TiT = (a+b,)-‘ |^2b2 “5^ +p + 4bj +a^4-^ + 5aj ■ Po 1,1+ da jj 
de ~ +l*l) |^^ 1*2 de de + de / Po de J 

^Pn r ^1*2 / \ 

= (a+b,)-‘ [2bj -5s + (4a-Po) as J ' 
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This may be written in the generalized form, for i = 1, 2, 3, 



From Equation 21.2, 


Thus, for i = 1, 2, 3, 


For convenience in later equations, we now differentiate the inverse of a„ p, = 
to find, for i = 1, 2, 3, 




da^ 


2«2 


(^oPo)'* = 


Po"' (a + bj)-* j^(a +bi)-i +Pg-1 ^ 


From Equations 31.1 and 31.2, neglecting terms of fourth order, 


u 


f2sy ( c2 \ 

y (i-s) sj 

(1-2S)J 


Now, for i = 1, 2 only, 


^ = - U»(l-S)|c» 4, (-oP.)"' - (I;) [f„ ("-ifF, ®)’5^ jI; (“iPo)’' 


Vf. 


'2S' 


<'-“> ■ ».P. 


(' ■ ^0 ") . 




d 

aq. 


(a+bj)po 
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However, 




a-5-o" 


- u*|c* [(1-S)^ (aoP„)-i-(a„p„)-»]-(|^) {[l“ifr„s][l 

+ c* (1 -S) (S ^ (aoPo)-‘ + (a, p„)-‘)|j^l + ^ (1 - 2S)j 

<1-S) (l--if^ s)[l+^ (1-2S)] [(1-2S)^ (a„p,)-*-2(a„p„)-‘ 


By Equation 32.1, 


^0 Po “ ■ 


Thus, for i = 1, 2, 3, 


fSi - 2 -i- I \-i + —sL 

aqj “ 5qj (^oPo) a^p, aqj 


Combinir^ Equations 32.1 and 32.2 gives 


Cl = £u(l-Cj(l-C,S)-‘ 


Po 


Thus, for i = 1, 2, 3, 


§ = I; (•-':.s)-'{(i-c,)(^ -^|^)->-(i-c,s)-[<i-S)^ -s„ (i-c,)c,]j 


Combining Equations 32.1 and 32.3 gives 


P = — u(l-S) (l-C,S)-> 


Thus, for i = 1, 2, 3, 


^ = p[( 1 -C*S)-»( 83 . C, + S^)-p-‘^ - 83 , (l-C,S)-‘^u 



Equation 21.1 may be rewritten as 


Thus, for i = 1, 2, 3, 


a, — “ 


5 /i(a +b,)' 



For convenience in later 

equations. 

, we now evaluate, for 

i = 1, 2, 3, 


A. lb. 


ab. 


\b. 

I = 

2 . 

d<ii 


a /ba 

d<ii \ P 

j = P-‘^-baP-^ 

dp 

aq. • 

Equation 5.31* gives 




A, = 

(l_e2)i/2 

n = 2 



Here, and in the following, P„ (x) is the Legendre polynomial with argument x of degree n. The 
definition of R is given by Equation 5.28, viz., r^ (x) = x" P„ (1/x), where 0 < x < 1. We shall de- 
note by P„' (X) the derivative of the Legendre polynomial with respect to the argument. Then, for 
i = 1, 2, 3, 


dAj 

5q7 


- Aj p"* 


Tq. ~ ^21 A, e(l-e2)-i+ p(l 


. 2 \ 1/2 


fe (t) £ "(t) Pn (^) R„-2 [(l - 






8,,pe 




♦Until otherwise indicated, equation numbers used for the defining relations will refer now to Reference 3. 
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By Equation 5.36, 




so that, for i = 1, 2, 3, 


(?Aj 


-»,p- ^ *2 p-‘4 |- (t) [(1-*")“] 

L n = l 




+ 8,, p-‘ e 


L n=l 


By Equations 5.61, 5.50, and 5.53, 


00 

A 3 = (i-e^)>/^p-3 2] D„R„,3 [(l-e3)V2] 


where is computed as follows; 




j 

= z; 


(_l)j-m (£> I 


pj 


(n an even integer), and 


m= 0 


(-i)i-'" 1-1 ' 


^2m + l V 


(n an odd integer). 
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Then, for i = 1, 2, 3, 




L n=0 

- £](n+2) [(l-e^)‘/*]"‘D„P„,, [(l-e*)' 


■ 1/2 


where d'D^/dq^ is computed as follows. If n is an even integer, then 


dH„ dD 


■3^ Sqj 




4 (T)i 


+ 2^l^i; (-l)'-"m(fl 1^1 


2m“l 


^ C \ 2 ( j "»")/ t >2 
VP j VP^ ^2mVb 


*4; z: © 


2(j-m) (h^ ^ I '-’I 


P2'4^j 


If n is an odd integer, then 


^ ^ '^Pai+i 

aqj aqj 


2cp-^^ ^(-l)i-"-(j -m) 

m= 0 


• Vp/ *’2ni+lVt>2y 


VtJ ? ,(-lV-"(2>n + l) ip] Vt] P2„«te 


3 (^ 2 ) ^ 


m= 0 


By Equation 5.37, 


^21 = (l -e2)i/2p-2 ejVjj + (3bj2 -b/)p ‘ - •§ b, b/p'* (l + ■? ®^) + "I ^> 2 '* P ^ (4 + 3e*)] 
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Thus, for i = 1, 2, 3, 
dk 


<jp 

dqj aqj 


- 2p-2 (3b/ -b/) ^ +2p-> (sbj ^ -bj 
+ I bj p-2 (l + e*) ^3bj bj P‘‘ - b^ ^ - 2b, 5^) 

+ -j bj3p-3 (4+ 3e2)^^ -bjp"‘ ^ 8j, b/ p'* e(b/p-‘ -b,)j 

By Equation 5,38, 

A 22 = i (l-e^)‘/^p-^e^ [(3b,^-b/) - 9b,b,^p-« +|b3^p-=>(6 + e^)] • 

Thus, for i = 1, 2, 3, 

-i r 

^2) 


aAj 

~dq 


)-3e2 

( 1 - 

y/2 ^ 


-> (3b,- 

5p 

+ ( 3 b, 

9 

-1 , . 


aq. 

- 2 P 

^2 

) dqj 



n 

1 ' 

1 

- 9p' 

•«b,b. 

-•^ 2 ] 


} 

e(l - 

e2)i/2 

[ 3 b,* 

1 

cr 

-9p-» 


By Equation 5.62, 

^31 ^ (1 “ P'^ e [2 + bj p‘^ (s+ e^) - p -2 ^2 + 0 ^) ( 4 + 3e2)j 


15 


Thus, for i = 1, 2, 3, 


aA 


■31 




3A3iP'‘^ + (l-e2)»/2p-4e||^-bjp-»(3+ |e*) 


+ 2p-* (4 + 3e2){4b/ + c2)j ^ + (a+ | e*) ^ 

- b^p-* (4 + 3e2) +S^. [| bj e-6p-‘e(| b/ + c2)]| 

*2i P’^ [(l e2 - (l -e2)i/2] j^2 +bj p"‘ (s + ^ e^) - p'^ b/ + c^) (4 + 3e2)J 


By the corrected version of Equation 5.63, 


A 32 = (l-e^)‘^^p'^e2[i+ Jbip-‘-p-2(|b/+c2)(| + J-e^) • 


Thus, for i = 1, 2, 3, 
dK 


'32 




- - A32 3p-‘ ^ + Sjj e(l -e 2)-i^ 


+ (1 - e^) >/2 p 4 g2 I p 2 1,^2 + c2 )( 3 + -^ ) 

3 , 1 ap 3 1 ^ 1 

- 4 P bjJ + 4 ^q. 2 P *^2 2 ^ y aqij 

+ Sjj (l-e2)i/2p-3ej| (1 +3bj p'*) -P'^(^ b/ + c2) (3 + e2)] 

By the corrected yersion of Equation 5.64, 

A 33 ^ 4 (1 ■®^)*^^P"‘ b, -p-i (I bj2 +c2)j . 

Thus, for i = 1, 2, 3, 


3qj 


- " 3 A,, p 


33P-‘^, (l-e^)‘/*p-^e3|[2p-^(jb/4c^) -^p-‘bj^ 


1 *^^2! 1 r / \ 

4 ^ " + 3 [3(1 -e") 


1/2 
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By Equation 5.65 


A34 = + 


Thus, for i = 1, 2, 3, 




*1*.. •’p-‘ [j (1-.’)-'" 

By Equation 47,* Q = ^ + s • Thus, for i = 1, 2, 3, 


_dQ 

dq 




From Equation 102.1, ^ - P(1 ~S) ^ , so that, for i = 1, 2, 3, 


^ = (l-S)-‘[^ - SaiPO-S)-', 


From Equations 147.1 and 147.2, 

hj " T ”^ 2 ) and h 2 - ~2 (l “^1 “^ 2 ) 

Thus, for i = 1, 2, 3, 



and 

^ = 2 h^f— + — • 

5qi ^2 V^Qi ^qj/ 

From Equations 100.1 and 100.2, 02 = Q(1~P)“^ and 63 = Q(l+P)”^ . 


♦Until otherwise indicated, equation numbers will now refer to Reference 7. 
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Thus, for i = 1, 2, 3, 


, ( dO dV\ 


and 




From Equation 131.2, e' = ae(a +b,)‘‘ , so that, for i = 1, 2, 3, 
ll: = (a+b,)-* j^Sj. e + Sjj a-e' ^8ij 


From Equation 149.1, (i-Su“*)*/^ , where sgna^ = agy/|aj| is +1 for a direct 

orbit and -1 for a retrograde orbit. Then, for i = 1, 2, 3, 


do... 




S [(l-Su-»)*/2 -S 3 .) 


By Equation 76, neglecting terms of fourth order. 


B,' = i Q* + - I C, PQ=* + I Cj P2 q 2 + ^ (4Cj + 3 Cj2) Q-* 


tIcjC^PQ^ + ^ (X3^ + 5C,*C3)q« + ^C/Q* 


Thus, for i = 1, 2, 3, 

+ +^C3Q^(X3+sc3*)]^ + 2P^ -|q(pq^ 

"2C.P^) +|pq(pQ^ +2C3Q^^ +2C3P^) + (2 3C, 

45 , / ^^2 dV dQ \ 

+CjPQ^ +C,C 3 Q 5 ^ +4C,C3P5^j 

1C r, dC^ 17*\ / dC,y <^50 \ 

+ ^ Q® |_{4C2 + 5C,2) ^ + IOC1C3 ^ (3Q^ +8C3 ^) 
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By Equation 65, neglecting terms of fifth order, 

B, = 1 -4ciP + + + -|C,C,PQ2 


-41 C,3PQ* + ^ ^ C*Q* 

-|c/P*Q^ -Tic.C/PQ^ 

Thus, for i = 1, 2, 3, 

^ +Ci ^.) + (I c," + 5 Cj) ( 2 P + <3 + (P" + I c, + 4 a^f) 

^Ac,Q3[..¥c.Q^^4(frc,^Q^](Q§^^a^)-|Q(c.-3^ 

+ CiPQ^ + 2C,C,P^) -t|c.^q(3PQ^ +2C,P^) 

+ mC,Q"(2C,Q^ +C,Q^ +4C,C,^) +4i C,3 q3(q^ + C, | C, PQ (pQ ^ 

+ +C,P^^) -41c,Q3(c,PQ^ + 2C,pq^ 

+ C,C,Q^ +4C,C,P^) + III C, Cj QS (c, Q + 3C, C, 

By the corrected version of Equation 95, neglecting terms of fifth order, 

Ba = -4c, -|c,3 - (^C,^+4fc,3c, + fc/)(l+4Q0 -|C/P^ 

- c,3c/ + ^c/)(i + 4q3 + |q^) - ^c/(i + 4q3 + |q^ + ^q«) 

+ Ci^"|c, c,)p + 4|c, c/(p+|pq3) 
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Thus, for i = 1, 2, 3, 


= -i(^*|c.^)-(.*4Q=)[|c.Sc>,x,)|^,|(|c.^.c.)^] 


r 35 


15 , 


3 

ci- 

45 


35 . 

( 

3 0 

15 .\1 

dQ 

3 

/ ac. 

“ [l28 

C ^ + 

16 Cl" 

C 2 

+ -8 

16 

C C 2 
'-1 ^-2 

P + T28 C/ 


■2 + 



4 C 2 P 


15 

/ 


3 


r? 



l^ , \ 

dC^ 

1 5 

/ 

3 


dC. 

- T6 C; 


2 + 

8 

Q^) 

[- 2 C, 

C 2 


bci"+cj 

(9q, 

J + 16 

C 211 + 

tQ") 

3P(^C, 


+ 2C^ 

dCA 




i 

f21 , 

\ 3q1 

35 

, / 

1 , 

3 ^ 

5 


d<ij 

+ 3Cj 

Cl 

! dq 

- - c 

i 

2 QI 

It Cl" 

^^2) 3qJ 

32 

C/(l 

+ 2 

+ 8 Q" 

^leQ' 

) <?qi 




By Equation 116.3, 


’®12 “ " 4 " 8 ^2 ’ 


SO that, for i = 1, 2, 3, 


dB, 


dq^ - 2 VVi-'-2 V ) ^ 

By the corrected version of Equations 116.4, 




B, 


21 


C 2 PQ + ^ C, C, Q3 + ^ Cj Q 


and 


B 22 " [(4C2+3C,=')Q2 + 3C,2Q^] 


Thus, for i = 1, 2, 3, 
aB, 


<^®2i 1/9 \ /9 \ ^^2 /27 „ 1 \ aO dB 

= yQli + lc^Q^)^ + q(^c,q^-p)^ + (^c,c,q^-C 3P+ -c^Q^ . 
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II 


and 


^^22 

d<i, 




+ Q(2+3C,Q2) ^ 


+ 3C.Q 


By Equations 122.1 and 122.2, 


2t7Vi = [-2cl^Y^^ [a +h^ + \ + Pi^^i » 

and 

2777^2 " ^2 ^2 ^2 * 


Thus, for i = 1, 2, 3, 


dq. 


{2nvi) = - 2t7v, |(a+b, +Aj +c2 AjBj'Bj i)-‘ |^S,j + 5 ^ 


3Pl. / 3A, aBj' aBjXl 1 

^ B,-‘ \b' ^ ^ -A, B/ B,-‘ ^ jj - 2 


and 


dq. 


( 2rrv^ ^ - 2777^2 U ^ ^ ^2 


[Sii 


ab, aA, 

+ a^ + a^ tc^Bj 


1 / . 

2 aqj 


aBj' aB 

\ aqj ~ ^2 ®2 * aq 


:)]-(' 


da2 _ ^^2 1 du 

2 ~ (9q. '*’ ^2 ^q. ~ 2 ^ ^q^ 


dA 


The following time -independent partial derivatives are used only when the differential 
rection includes periodic terms throu^ the second order. 

By Equations 5.32, 5.33, 5.39, and 5.40 of Reference 3, 

All = ^ ( 1 -e2)i/2p-3b/ e(b/ -2b, p) , 


'12 “ 32 


(1 -e2)»/2p-3b4 e2 


A 23 = i (l-e^)‘/^p'"bj 2 e3(b2*p-‘ -b,) 


aa, 


sq, 


cor 


and 


'24 " 256 


(l_e2)i/2p-Sb4e4 . 


Thus, for i = 1, 2, 3, 


3A 


11 




- Aj, jsp-i ^ +S 2 i [e(l-e2)-» -e-»]| 


(l-e2)i/2p-3bje[b, bj ^ +^2?^ +2(bjp-b/)^j 

^ = I (1 -e3)>/2 p-3 b^3 e3 (^ - I p- S 3 , p-3 eb 3 ^ [(l -e3)V2 _ ^ (1 -e3)-i/2 e*] 


aA 


7 = - A23 P'* :^ + I (l-e2)i/2p-‘*b3 e3|^b3p-i{3b, -4b33p-») 

abj ^^2! 

- t>2 ^ +2(2b33p-l-b,) 


and 


dA. 


+ ^ 82 , p '•bj^e^ (bj^p-i-bj) [3(l-e3)»/2- (l - e^] 


7 - Aj4 (4b3-‘ -5p-‘ ^ 83, p-Sb3‘»e3 [4(l-e2)i/2_ (i_g2)-V2e2] 


dq 


By the corrected versions of Equations 116.3 and 116.4 of Reference 7, 


B„ = - 2PQ + 8 


®13 ~ 24 ’ 


Bi4 = mCjQ" 


®23 “ " 16 ^2 ’ 
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and 


®24 “ 256 ^ 2 ^ • 


Thus, for i = 1, 2, 3, 


and 


aB„ 

(9 , 


" U^iQ^ 

<?Bi3 

1 , 

aq. 

= -mQ" 

aBj, 

1 , / 

aq; 

= 

^®23 

1 , 

dq. 





(9 , \ dO aP 3 , 

(8C,Q*-2p)^^ - 2Q^ ■ 

1 ,( 3Q\ 

- +3C, ■ 

1 ,( ^^2 3Q\ 

1 / ^Cj aCj ^Q \ 

+3CiC, 


128 ^2 


J ag\ 


TIME-VARYING PARTIAL DERIVATIVES IN THE DIFFERENTIAL CORRECTION 


We shall continue the use of the generalized notation introduced in the preceding section, 
whereby we let qj(i = 1, 2, 3) refer to the orbital elements a, e, and S, respectively. The 
time-dependent parameters will involve, additionally, the orbital elements /3j, and /3j, which we 
shall represent as /3j (i = 1, 2, 3) in the generalized form of the partial derivatives. 

By combining Equations 121.1* and 123.1, 

” 2ttv, (t +/3j - c^/3j aj”! Bj' Bj"^) . 

Thus, for i = 1, 2, 3, 

aM^ , 

Wi " (t (2 t7Vi) 


- 277T/j C* /8j ttj * Bj * 


/as,' aa^ 

V'a^ -®i'“2'* 




a^\ 

d9i) 


•AH equation numbers used in specifying the defining relationship for a given variable will henceforth refer to 'Reference 7. 
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Also, 


= 27n/, . 


dill^ 

W2 " - 27TV, 


and 


<?Ms 

W, = 0 ■ 


In what follows, whenever a partial derivative with respect to /3^ is zero, it will not be given. 
By combining Equations 121.2 and 123.2, 


'•'s - 277^2 [t +^J V* . 


Thus, for i = 1, 2, 3, 

- [t +;S, +/?2 Aj-I (a + bj + Aj)J (277^2) 


Also, 


+ 2rrv^P^a-^ 



abj 


i?Aj / da^ dk^ 

+ ^ - (a+bj +Aj) ^a2-» ^ +A2-« ^ 


^'J's 

^ - 277^2 , 


and 


d;S 2 " 2777^2 a 2 -» A 2 -* (a + bj + Aj) . 


By Equation 131.1, letting 8 s +E^ gives 


8 - e ' sin 8 - M . 

s 
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Thus, for i = 1, 2, 3, 


dS 


( 1 - e ' cos S) * 




+ 


dcij ‘ 


Also, for i = 1, 2, 


d8 _ , D 1 

d/3. - ( 1 -^ cosEy^ ^ . 


By the anomaly connections given in Equations 132, taken to zeroth order, 


cos v' - (cos S ** e) ( 1 - e cos 8) ^ , 


and 


sinv' - (l - e^) ( 1 - e cos 8)“^ sin 8 , 


where 


V 


0 


v' “ M . 

s 


Then, for i = 1, 2, 3, 

^Vq r aS 1 

^ = Kl - e^) sin S ^ + §21 sin^ Sj (sin v')“‘ (1 - e cos S)'2 - . 

Also, for i = 1, 2, 

^ = (l-e = ) sin8(sinv')-‘ (l-ecosS)-2 ^ ^ . 

By Equation 133, 

'^0 = ^2 . 

Thus, for i = 1, 2, 3, 


dq, 





da- 2 

Vo ^ - 


1 -1. *^''1 

2 “2 “i Aj Vo 


1 -1 . du 

2 “2 “ Vo - 5 ^ 


dP^2 

^2 dq^ " ^2 ^2 ®2 ^ aq, *** ^2 ^2 ~d^^ 


t) 
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Also, for i = 1, 2, 




1/ * n V 

W[ " “2 ^ • 

By Equation 136.2, 

Ml = - (a +b,)-> [(Aj +c 2 A 2 Bj' B^- i) Vo + c^a/* (- 2 ai)*/ 2 ul/ 2 B^^ sin 2 {¥, +S»o)] . 

Thus, for i = 1, 2, 3, 

aM f / abj\ av. faA, / aA, 

^ - - (a +bi)-‘ |m, ^§11 + + (Aj +c2 AjBj'Bj-i) - 5 ;^ +Vo[^^ +c2Bj-> ^Bj' 

<9B * VI /(?B 

-AjBi'Bj-i ^jj + c^a/i (- 2 a,)i/ 2 u»/ 2 sin 2 {¥^+S'o)(^- 5 ^ + I ^ 

1 ^“1 ^“ 2 \ 

- ^ “f* Bi, ^ Bi, . 2 c^ a,-‘ (- 2 ai )>/^ u*/^ Bi, cos 2 («P, ^ 

Also, for i = 1, 2, 


dM 


1 r 

- = - (a+bi)-» j^(Ai +C 2 A 2 B 1 ' Bj-*) ^ + 2 c 2 a/» (- 2 aj) »/2 ui /2 cos 2(V^ +^o)(^ * ^j\ 


By Equation 137, neglecting terms of third order, 


Ej “ ( 1 - e ' cos £)“ ^Mj -^e'(l~e' cos S)“ ^ sin 8 


Thus, 

for i = 

1, 2, 3, 










aEj 




' COS 8)“ ^ 

( 

a£ 


de'\ 



' cos £)”® 

^Qi 

(l-e' 

o 

* a^ " 

Ml (l-e 

le ' sin c 


- COS 8 

5qJ 

-M, 

(l-e 



3 


o 

0 

U) 

00 

1 

/, • o 

a£ 


ae'\ 

1 


o SS 



- 2- Ml 

e' (1 -e' 

fe sin 8 

aqj 

- COS 8 

aqj 

2 

Ml e' 



1 ^ 
2 “i aq, 


Also, for i = 1, 2, 


aE aM 

= (1 - e' cos Sy ^ [l -Mj e' (1 - e' cos S)"2 sin Sj 


^ g r 1 3 “] 

Mj e' ( 1 - e ' cos £)" * |sinS + ’^Mj (1-e' cos S)” * cos S-" 2 Mie'(l-e' cos £)” * sin^ Sj 
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By the anomaly connections, taken to first order, 


and 


where 


cos v" - jcos (S+Ej) ““ ej - e cos (S+Ej)J ^ , 


sinv" - [l - e cos (S + E jJ ^ sin (£ +E^) 


Vi - v" - v' - v" - (M^ + Vq) 


Then, for i = 1, 2, 3, 

(9q. 


(l-e2) sin (e+Ej) 


+ §2i sin^ (£ +Ej) 


(sinv")-i j"i - ecos (S+E,)] * 


Also, for i = 1, 2, 


= (l -e2) sin (8+Ej) (sinv")-i [l - e cos (8 +Ej)J ^ • 


In the foregoing and in the following, we employ the symbolic notation 


and 


(9 V ' 

5q. 


dw' 

w: 


5Qi 


a/3. + a/3, 


By the corrected version of Equation 139, 


'Pi = -Bj * [b ,2 sin 2{'P^ cos (V, +¥,)] 


+ aj(-2aj) Bj" ‘ (Aj Vj + Ajj sin v' + Ajj sin 2v') 
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Thus, for i = 1, 2, 3, 


d'V dB. f dB 

= B;* [b,, sin 2 («P, +¥,) + B,1 cos («P, +<?jj ~ B; ^ jsin 2 («p, +>pj 

5Bjj r T /^'•'s 

tcos («P^ +¥„) ^ + [2B,, cos 2('P, +«P„) sin («P^ +»P,)] + j^jj 

+ (-2aj)- »/2 u- B; » ^ - 2 “2 (“f * a^ + " a^ + 2 B - 1 (a^ Vj 


+ Aj j sin v' + Ajj sin 2v') + 


r ^ 

['i 


avj 

+ +sinv'-^ 


V ^^2 2 5 V ^ I 

+ sin 2 v' + (Ajjcosv' +2A2 jCos2v') 

Also, for i = 1, 2, 

5*Pj P /av a>Po\ 

W, " - ^2 ‘ P ^22 H'^S + 'Po) - J Sin (<P^ + «P„)J 

[ <9 V - n /“I 

^2'^+ (AjjCosv' + 2 AjjCos 2 v') 

The following time-dependent partial derivatives are used only when the differential correc 
tion includes periodic terms through the second order. 

By Equation 143.2, 

M 2 = -(a+bj)”*^ ‘^AjVj + Ajjsinv'+Aj2sin2v' 

+ C*a;i (-2aj)i/2u»/2 [b^> +BjjCos (v, + <Po) 

+ 2B,,<PjCos2(>P^+«Po) +B,3Cos 3(«P^+<P,) + B,, sin 4 («P^ + «P„)] j. 
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Thus, for i = 1, 2, 3, 

<9M 

= ' (a + b^)-» fs, 


ab; 

\ 

aAj 


avj 


aA,3 


aA,3 

cr 

4. 

J.V 

1 aq. 

Ai 

aq, 

sin V ' 

aq, + 

sin 2v' 

aq. 

2A,3 

cos2v) ^ 

+ c 


(-2“i) 

1/2 1/2 


au 

aq. 


aa2\ 

r 







“2 * 

aqij 

[b/ -p. 

+ B 

11 COS 


0 ) +2B 

12^1 cos 

2('P, +«Po) 


as; 


''i aQi 


+ B,3Cos3(vP^+¥„) +B,^sin4(V^ + V„)] + a~ ^ 

+ Bi' ^ K+n) ^ + -^ + ®,2^j cos2(<l«^+«pJ 

aBj3 aB 

+ cos3(l»^ +V„) +Sin4(s»^+«P„) - 3 ^ - [Bjj sin +4«o) 

+ 4 B 3 3 <P, sin 2(<P^ +<!<,) + 3 B ,3 sin 3(<P, + 'P„) - 4 B,, cos 4(«|.^ + -P,)] 


Also, for i = 1, 2, 


aM, 


T dv 


dv' 


dfi. - -(a+bj)-MAj^ + (AjjCosv' +2 Aj3Cos2v') 


+ c* ttj * (-2aj) 


a«p, a«p 

W, +2B,3Cos 2(>P^+V,) ^ - [Bjjsin (>P^ + «pj 


+ 4B33^-3sin2(<P^+«P„) +3B,3sin3(<P,+<P„) - 4 B,, cos 4(¥^ + <pj] ^ J 

By Equation 143.1, 

Ej = - e' cos (£ + E,)J . 

Thus, for i = 1, 2, 3, 
aq, = [l-e' cos (8 +Ej)J ^ 


- M, 


[l-e'cos(8 + E,)]-^ e'sin(8 + E3) ' cos (8 . E J 
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Also, for i = 1, 2, 


dP; [l e'cos(S+Ej)] Mj [l -e' cos (8 + Ej] sin (e + Ej) ^ 


Since E = S +E, tE^, then, for i = 1, 2, 3, 


Also, for i = 1, 2, 


_£E _ ^ ^ 

<9qi 


aE _ ^ 

dp. - dp. + a^S, + dp. • 


By the anomaly connections, taken to second order (where we choose the notation v rather 
than V "'), 


COS V - (cosE-e) (l-ecosE)“ 


and 


where 


sinv - (l - (l - e cos sinE 


- V - V = V - (M +v 




Then, for i = 1, 2, 3, 


- j^(l -e2) sinE ^ + S^. sin^sj (sinv)-i (1 -ecosE)-2 - 


Also, for i = 1, 2, 


avj 

W 


- (l -e2) sinE(sin v)"»(l-ecosE)-2 ^ 


Since v = v' +v, +Vj , then, for i = 1, 2, 3, 


Jz = ^}l 1 fh 

^<li ~ dqj ^ aqj aqj 
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Also, for i = 1, 2, 


dv ^ ^''i ^^2 

W[ + d/3. + dfi. 


By the corrected version of Equation 145, 


V, = 


[■Bji'J'iSin + +2B,j'P,cos2(«P^+¥o) + cos 3 (<P^ + «pj 

+ Bj^sin4(s«^ +Wo)] +aj {-2a^y u~ ^ (A^ j Vj cos v' 


+ 2 Aj 2 Vj cos 2v' + Aj 3 sin 3v' + A^^ sin 4v') 


Thus, for i = 1, 2, 3, 


lif = [-B3^>PjSin(«P^+«Po) +2B„«PjCos2(«P,+'Po) +B,3 Cos3('P^ + 'pJ 

as f / 5B a^,\ 

+ B,, sin 4(«P^ +«?„)] +B 3 ,^j sin(<P^+«pJ 


+ 2 


/ aBj2 a?>\ 

(;Pi ^ cos2('P^ 


as. 


aB 


+ V„) +cos3(¥^ + Vo) +sin4(^^ +?>o) 


[Bj , cos +?<„)+ 4 B32 sin 2(^3 + +3 B^j sin 3(^-3 + ¥„) 


- 4Bj^ cos 


aB 

^2b;* 






1 / _ , _ - du 


(^2 ^2 ^ ^2 1 ^ ^ ^ ^22 1 ^ ^23 sin 3v ' + Aj 4 sin 4v 

+ a, (-2a,)->/*u->/2B;i |v, ^ + A, ^ + A,i 

/ aA^^ avj\ aA aA 

•"2^v, cos2v' +Sin3v' +sin4v' 


av' 


(Ajj Vj sin v' +4 Ajj Vj sin 2v' - 3 Ajg cos 3v' “ 4 A^^ cos 4v') 
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Also, for i = 1, 2, 




1 +’*' 0 ) - 2 B 22 COS 2(w^ +W„)] ^ + [Bjjfj cos +9<o) 

/d'^ 

+ sin2(-P, + -P„) +3B,3sin3(?r^+-|-^) - 4 B,, cos 4(¥^ + ¥^)] 

P ^^2 

+ a 2 {- 2 aj)-i/ 2 u-V 2 Bp _ + (a^ ^ cos v ' + 2 cos 2v') ^ 



(A 2 iViSinv +4A2jVjSin2v' 3A23Cos3v'~4A2^ cos 4v ') 

Since W = + 'P(, + "Pj + Wj , then for i = 1, 2, 3, 


^ ^ ^ ^ ^ 

<?qi i3qi i5qj ■*■ 3q. 3qj 


Also, for i = 1, 2, 


^ ^ ^ ^ 

d/3. - 3/3. + 3/3. + 3/3. + 3/3. 


This concludes the evaluation of the partial derivatives of the uniformizing variables E, v , 
and 4* when the computation is followed through terms of the second order. If, however, second- 
order precision is not required, then the partial derivatives of , E^, v^, and may be omitted, 
and the foregoing partial derivatives of the uniformizing variables reduce to 


3E _ 

' 3qj ^ <?qi 


3v ^ ^ 

3qj (9qj 3qj 


and 


3^ 3'V^ aw, 

5qj 5q. 3q. (Jq^ 


The partial derivatives of E, v, and 4< with respect to (3.^ (i = 1, 2) are similarly reduced. 
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We now continue with the necessary equations preparatory to determining the partial deriva- 
tives of the inertial coordinates X, Y, andz. 

By special cases of Equations 104 with y = 'l'±(l/ 2 ) 7 r, 


and 


cos E2' - (1 - 02 sin 'p)" ^ (0 2”sinW) , 

sinE2 = (1 - e 2 sin v) “ ^ cos ^ , 


cos E3 " (l+e^sin^)“^ (e^+sin^) , 

sinE3 = (1 + 03 sin w)“ ^ cos ^ . 


Thus, for i = 1, 2, 3, 




Ti; - (l-e2sinS>)-i 


, , 1/2 as* , 


5 ^ 2 ! 


and 


3E ' 

~a^. = (l + ejSinV)'! 


. O \ 1/2 d'V 


5q. " 


V2 ^ 

dq, 


Also, for i = 1, 2, 


5Ej' 

~dW, 


/ \ - 1 / o\ 1''2 I? S’ 

( 1-63 Sins') (l-ej 2 ) 


and 




/ ^-l , ,vl/2 5S< 

( 1 + 83 Sins') 1 ( 1 - 832 ) ^ 


By Equations 114.1 and 114.2, 

Yo = i (1-20'‘^"E3* + I (1 + 20‘‘^^E3' 


and 

Xi = J (1-20‘‘''"E3' - I (1 + 20‘‘^'E3' • 
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Thus, for i = 1, 2, 3, 

aXr, 1 5E,' 


oxn 1 '"='2 1 1 r -> A r 

^ = 2 ( 1 - 20 '*"=*^. + 2 + +2 [ e 2 '( 1 - 20 -^=* " (1 + . 


and 


5qj 


Kfl 1 ^E' , dJ£- , r ,, 

- 2 (1-20‘*^*-5^ - 2 a + 20‘*^*-5^ +1 [e2'(1-2C)‘^" + E3'(1+25)'^^*]-^ • 

Also, for i = 1, 2, 


aE,' 


dE^ 


^Xn 1 aEj* » 

Wi " 2 <1 ‘20‘*'* -^ + J (1 + ’ 


and 


^Ej 1 


^Ej' 


a/3j " 2 20 ays. 2 (i‘‘’20 a/3. 

By Equation 150, 

4> = ySj - c^ttj (-2aj)‘*/^ (A3v+A3,sinv + A32sin2v+A33sin3v + A3^sin4v) 

+ tt3 a- » |( 1 - S)- 1/2 |^(hj + hj) Xo + (hi - hj) xj 


3 3 1 

+ BjV- 4 CjC^QcosV + 32 C32Q2sin2«'| 


Thus, for i = 1, 2, 3, 


^ = -c2 |^(-2a,)-i/2^ +a3(-2aj)-3/2-5^j (A 

^ r ^^3 

, 4 v) - c2a3 (-2 aj)-i/2 jv ^ 


3 V + A3 j sin V + A32 sin 2 v 


+ A33 sin 3 v + A3^ sin - 


dAai *^^3 2 

'^sin2v-^^ 


<9A33 ^^34 

+ sin 3 v — + sin 4 v + (A3 + A3 j cos v + 2 A32 cos 2 v + 3 A33 cos 3 v 


+ 4A3^cos4v) ^ 


N) Xo 


““3“2‘‘ aqj +T“3l^"* [(hj + 

3 3 ''i 

+ (hj - hj) Xi] + Bj ¥ - 4 C, Cj Q cos S? + 32 sin 2¥| 

+ a 3 a^- 1 ui/=i || S 3 . (1 - S)-^/^ [(hi +h3) Xo + (hj - h^) xj + d " S)' »/* |jhj + h^) 

axi 5 hj ahj] aBj 3 / aCj 

+ (hi-hj) ^ + (xo + Xi) ■^ + (xo - Xi) "^J aq, " T ^ 

^^2 aQ\ 3 f ^^2 aQ\ 

^CiQ^ +CiC 3 ^j cos«P +36 C 3 Q(q^+ sin 24- 

/ 3 3 , , \ a^-l 

+ (^3 + 4CiC,Qsin4- + 16 Cj^Q^cosOTj-^^ 


34 



Also, for i = 1, 2, 


~ (^3 ^31 2 A32 cos 2 v + 3 A3 3 cos 3 v + 4 cos 4 v) 


r r ^^0 

|(l-S)-»/2 [{h,+hj) ^ + (h,-hj) 


( 3 3 \ 

®3 + + Y6-Cj2Q2cos2<pj ^ 


Since the spheroidal coordinates are given by 


p - a (1 - e cos E) 


77 - P + Q sin^ 


1 - e cos E - a s 


^e\ 


Also, for i = 1, 2, 


a [ e sin E ^ - cos El 


Js ■ a e sin £5^ 


- a e sin £3^ 
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Further, for i = 1, 2, 3, 


dv dP . dO d'S 

^ ^ + sin«P^ + Pcos-P^ . 


and, for i = 1, 2, 


dp. dp, 


The partial derivatives of p and v with respect to P^ are both zero. 
Finally, the inertial rectangular coordinates are given by 


X - ]/(p2 + c^) (l - T)^) cos 4 > , 
Y = '}/{p^ + c^) (l - 7)2) sin<^ , 


and 


Z - pT) - h . 


Thus, for i = 1, 2, 3, 


dX 

dq, 


X 


p(p^■^c^) ‘ 


lE 

<?qi 


77 (1 - V^) 


JjL 

dq.^ 


d4> 


dX 

dq, 


p(p^ 


+ c 


) ^q, ■ ^(1-^") ,9q. 


+ X 


(9<^ 


and 


dZ _ d-rj dp 

dq, dq, 


Also, for i = 1, 2, 


<9X 

dP, 


X 





dr) 

dP, 
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and 


and 



17 (l - 77* 


.-1 ^ 

) dP. 


d<f> 

dP, 


dZ _ dn dp 

dp. - P dp.^^ V dp. 



dY 

^ = X . 


dZ 

JT, " 0 


OBSERVATIONS OF THE RIGHT ASCENSION-DECLINATION TYPE 

The differential coefficients in the form of partial derivatives of the inertial rectangular 
coordinates with respect to the orbital elements are completely general in the sense that they are 
functions only of the mathematical theory of orbital satellite motion. Thus, they are applicable to 
observations of spacecraft position recorded in any format whatsoever. Previously, it was 
assumed that the observational data were recorded as direction cosines with respect to a topocentric 
latitude-longitude-zenith coordinate system. Another format frequently used for recording obser- 
vational data is the right ascension-declination type. In this section, we shall discuss the minor 
variations in the equations that arise when this type of data is utilized. 

The coordinate system adopted for the use of right ascension-declination data is also situated 
at the tracking station on the earth^s surface, but its three coordinate axes are parallel to the 
respective axes of the inertial system. Again designating the topocentric ’local^’ coordinates by 
the subscript in this case, the Zj^^-axis is parallel to the earth^s polar axis, and the \ 
plane is parallel to the equatorial plane of the earth. The X„ -axis extends toward the vernal 
equinox, with the Y„-axis orthogonally to the east to form a right-handed system. 

The observed right ascension % is measured eastward from the vernal equinox, and the ob- 
served declination is measured as positive above, and as negative below, the equatorial plane. 
The corresponding computed values of the right ascension and the declination are given in terms 
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of the local coordinates by 


arctan 


(»■ 


and 


8 


arctan 




2+Y 

M } 


1/2 


In order to obtain a satellite’s local coordinates from its inertial coordinates, the inertial 
coordinates of the observation point at the time of observation, denoted (X^, Y^, Z^), must be 
known. However, no rotations are necessary to bring the two systems into coincidence in this 
case, since the topocentric and inertial coordinate systems are parallel. A single translation will 
suffice. Hence, the relationships for the local coordinates of the satellite are simply 


X ♦ , 

- Y - Y^ , 


and 


Zm = Z - . 

Notice that the foregoing simplified relationships are obtained from those of the direction-cosine- 
data case by the artifice of setting = 0 and 0^ = tt/2 in the corresponding equations for x^, , 

Y*, , and given earlier. 

The first-order Taylor’s series expansion for the equations of condition corresponding to each 
time of observation are 


and 



i= 1 



i= 1 
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where q. (i = 1, 2, • • • , 6) are the orbital elements. Expanding the foregoing partial derivatives by 
the chain rule yields 


da da dX^. da dY^ da 

and 

^§c aXj, dS^ dY^ dB^ dZ^ 

(jQi “ d\ 5q. + 3q. + 3Z„ 5q. 


From the expressions for and in terms of the local coordinates, we find - 


d a 

C 

\ 



da 

C 

\ 


’ 

d a 



0 , 


x„z„ 

~ 

^mV 

dK 

Ym Zm 

■ 



IlL 


where R„ - and r„ - ( ^ + y^ 2 ) 

Since the station coordinates x^, Y^, and z^ are independent of the orbital elements (and 
merely geodesic functions), the following simple relationships apply. 


^ _ J]L 

i?q. 
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and 



dY 


aqj • 


az 

aqi 


j , and dz/aqj 


= 1, 2, • • • , 6) are precisely those 
that have been evaluated previously in the differential correction scheme. 


REMARKS 

The differential correction process removes inaccuracies of the initial conditions (the nominal 
observations) and accounts for the effects of forces not considered by the analytical orbital theory. 
Such neglected forces may include aerod 3 mamic drag, electromagnetic effects, solar radiation 
pressure, meteoric bombardment, and residual gravitational influences (including those arising 
from lack of spherical symmetry in the satellite, as well as perturbing planetary potentials). This 
is all accomplished by producing a mean set of orbital elements throu^ an iterated least- squares 
fitting of the first-order Taylor's series expansion of the conditional equations to numerous ob- 
servational values. Generally speaking, the fitting will improve as greater numbers of observations 
are considered and as the time span represented by the observational data is lengthened. However, 
the complexity of the mathematical processes involved in the fitting increases rapidly as additional 
observations are admitted. Because of this latter constraint, it is often advisable to perform re- 
peated differential corrections at various intervals of time (known as "epochs") rather than at- 
tempt to accommodate all of the data in a sii^le iterated fitting. 

The orbital improvement method of differential correction discussed in this paper has been 
programmed, primarily in the FORTRAN language, for use on a large-scale electronic digital 
computer. The analytical nature of the entire procedure assures a very rapid computational 
process. Application of the differential correction (combined with an orbit generator of position 
and velocity components) has been made to botii actual observational data of artificial earth sat- 
ellites and to artificially generated "data" for extremum cases, e.g., polar and equatorial orbits. 
The results have proven entirely favorable. Particular experimental applications will be published 
in a later paper. 
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Appendix A 

Functional Dependencies of the Partial Derivatives 

The following tables are intended to display the interrelationship of the various parameters 
whose partial derivatives appear in the differential correction. The partial derivative of each 
parameter in the left column is seen to be a function of those, and only those, partial derivatives 
of parameters occurring in the respective line of the right column. There is a certain amount of 
flexibility in the ordering of the partial derivatives occurring in the differential correction, as is 
demonstrated by the functional dependencies illustrated by these tables. 

Note: In Tables A1 and A2, the asterisk indicates partial derivatives of the parameters used 
only if the differential correction includes periodic terms through the second order. 

Table A1 

Time-independent Partial Derivatives 
(All taken with respect to orbital elements a, e, and s.) 



Functional dependence on other 

Partial derivative 

partial derivatives 

p 

(none) 


p 

b. 

P. b, 

Po 

P. bj, bj 

“2 

Po 

(^0 Po)‘ ‘ 

b,. Po 

u 

Po • (®o Po)' ‘ 

C2 

(^0 Po)'*- 

Cl 

Po- «- Cj 

P 

Po - Cj 

“1 


fe) 

b,- b^ 

(^) 

P. bj 

Ai 

- (S (;) 
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Table A1 (Continued) 


Partial derivative 

Functional dependence on other 
partial derivatives 


P. bj 


Cj, P, Q 


Cj, Q 


C,,Q 


C,.C,.Q 


C^, Q 


Table A2 

Time-varying Partial Derivatives (All taken with respect to orbital elements 
a, e, S, and Exception: 4 >^ X, and Y are taken with respect to /33, as well.) 


Partial derivative 

Functional dependence on other 
partial derivatives 

Ms 

a^, Bj', Bj, (2nv^) 

fs 

b,, aj, Aj, Aj, {2nv^) 

s 

e', 

''o 

M , 8 

s ’ 

'i'o 

“a- ‘^l’ ^2’ ® 2 ’ ''o 

Ml 

t>j, u, aj, Aj, Ajr Bj , B 2 » Bj 2 » '^o* ^0 ^ 

El 

e' , 8, Mj 

V ' 

Ms- Vo 

''l 

8, Ej, v' 

'Pi 

o-j, u, Aj, Ajj, Ajj, Bj, Bjj, Bjj, 4<(,, v , Vj 


bi, tij, u, ttj, Aj, Bj , Bjj, Ajj, Ajj, Bjj, Bjj, Bj^, 


> 

> 

O 

*Ej 

e', 8, Ej, Mj 

E 

S, Ej, Ej 


v\ Vj, E 

V 

v', Vj, Vj 
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Table A2 (Continued) 


Partial derivative 

Functional dependence on other 
partial derivatives 


^2* ^2’ ^21* ^22* ^2’ ®21* ®22’ ^23’ ^24* 



’P 

"Ps. '•'o. '?!• 'P2 

e; 

e2» W 

e; 

ej. ¥ 


E^' . E3' 


5. Ej'. E3' 


a^, u, Cj, ttj, A^» A 22 » A 33 , A 3 ^, Q, hj, 


h,, a3,B3, v,y,Xo* 

P 

E 

V 

P, Q, ^ 

X 

P^ V 

Y 

4>^ p^ V 

Z 

Py V 


46 


NASA -Langley, 1967 


30 



^'Tbe aeronautical and space activities of the United States shall he 
conducted so as to contribute to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— ^National Aeronautics and Space Act of 1958 
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